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ABSTRACT 

Aims. NGC 2903 is a nearby barred spiral with an active starburst in the center and H ii regions distributed along its bar. We aim to 
analyse the star formation properties in the bar region of NGC 2903 and study the links with the typical bar morphological features. 
Methods. A combination of space and ground-based data from the far-ultraviolet to the sub-millimeter spectral ranges is used to 
create a panchromatic view of the NGC 2903 bar. We produce two catalogues: one for the current star formation regions, as traced 
by the Ha compact emission, and a second one for the ultraviolet (UV) emitting knots, containing positions and luminosities. From 
them we have obtained ultraviolet colours, star formation rates, dust attenuation and Ha EWs, and their spatial distribution have been 
analysed. Stellar cluster ages have been estimated using stellar population synthesis models (Starburst99). 

Results. NGC 2903 is a complex galaxy, with a very different morphology on each spectral band. The CO (/=l-0) and the 3.6 fim 
emission trace each other in a clear barred structure, while the Ha leads both components and it has an s-shape distribution. The UV 
emission is patchy and does not resemble a bar. The UV emission is also characterised by a number of regions located symmetrically 
with respect to the galaxy center, almost perpendicular to the bar, in a spiral shape covering the inner ~ 2.5kpc. These regions do not 
show a significant Ha nor 24 fim emission. We have estimated ages for these regions ranging from 150 to 320 Myr, being older than 
the rest of the UV knots, which have ages lower than 10 Myr. The SFR calculated from the UV emission is ~0.4 M© yr~\ compatible 
with the SFR as derived from Ha calibrations (~1 M© yr~0- 

Key words. Galaxies: barred - Galaxies: individual (NGC 2903) - Galaxies: evolution - Galaxies: ISM - Galaxies: dynamics - ISM: 
HII regions - Ultraviolet: galaxies 



1. Introduction 

In the years after the launch of the Spitzer Space Telescope and 
the Galaxy Evolution Explorer (galex) knowledge on the star for- 
mation (SF) in galaxies have grown considerably. The panchro- 
matic view of nearby galaxies offered by large surveys carried 
out from these space telescope missions gives an extraordinary 
database to study star formation in galaxies. In particular, it al- 
lows us to link general galaxy properties with the the local inter- 
stellar medium (ISM) properties and the galaxy dynamics. These 
links are crucial for understanding SF in galaxies. In nearby 
galaxies these sets of multi- wavelength data give the opportunity 
to study in detail, with good spatial resolution, the location and 
properties of relatively young populations, recent massive star 
formation and dust attenuation (e.g. Calzetti et al. 2005; Tamura 
et al. 2009; Relano & Kennicutt 2009). 

New dust and SF indicators have been developed and 
calibrated using all the available bands. The combination of 
24 pm and Ha emission has become a reliable indicator of dust 
attenuation (Kennicutt et al. 2007; Relano & Kennicutt 2009). 
Both bands are linked to star formation processes, with the 
Ha emission originating from the recombination of hydrogen 
in the surrounding medium of very recently formed (less than 



Table 1. Log of ground-based Ha observations. 



Band 


Date 


Filter 


Exp. Time 


Seeing 


Ha 


29 Oct 2007 


WFCH6568 


3xl200s 


1.4" 


R 


30 Oct 2007 


HARRIS R 


3x250s 


2.1" 



a few Myr) massive stars, and the 24 pm emission as tracing 
local star formation radiation obscured by dust (Calzetti et al. 
2005). Extinction corrected UV emission can also be used to 
retrieve the star formation rate (SFR) (Kennicutt 1998) and to 
get some insight into the star formation occurred in last Gyr 
(e.g. Bianchi et al. 2005; Hibbard et al. 2005). This feature 
makes a combination of UV, Ha and 24 pm emission ideal 
to reconstruct the recent star formation history in a galaxy by 
studying the location and properties of individual clusters and 
those of the gas and dust emission around them. 

To this date, and to our knowledge, only a few galaxies have 
been analysed in this way; M51 (Calzetti et al. 2005), M81 
(Perez-Gonzalez et al. 2006); M33 (Relano & Kennicutt 2009; 
Verley et al. 2009); NGC 7331 (Thilker et al. 2007). 

Barred galaxies offer a useful tool to investigate the physical 
conditions that favour star formation in galaxies. The motions in 
bars are characterised by non-circular motions that push the gas 
into intersecting orbits where shocks and star formation can be 
triggered. The position and strength of the shocks are determined 
by the bar potential, the global dynamics within the bar region 
is driven by the bar. We have now a relatively good understand- 
ing of the gas behaviour under a bar potential (e.g. Perez et al. 
2004) and this knowledge can be used to understand the condi- 
tions triggering star formation. Bars and their surroundings host 
extreme physical conditions and a variety of ISM environments. 
They are perfect places to study the link between the conditions 
favouring star formation with the galaxy dynamics. Therefore, a 
panchromatic view revealing the history of star formation in bars 
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can give a unique insight into the links of star formation and the 
galaxy dynamics, which could also help to understand how bars 
form and evolve. 

This work presents a detailed multi- wavelength study of the 
star formation in the bar region of NGC 2903 analysing the cor- 
relations between the location and ages of the young stellar clus- 
ters with the morphology of the bar. This is done by analysing 
the emission in Ha , UV and 24 jim as well as optical data 
from the Sloan Digital Sky Survey (SDSS), and using 8 yum and 
CO (J=l-0) emission as complementary data. 

NGC 2903 is chosen for this research for a number of rea- 
sons: it is close by (8.9 Mpc; Drozdovsky & Karachentsev 2000) 
allowing us to have high spatial resolution (~ 43 pc arcsec"^), 
and it is isolated from large companions, preventing major 
merger effects in the results. Recent work (Irwin et al. 2009) 
on the H i content of NGC 2903 has shown that it possesses a 
large H i envelope of around three times its optical size. They 
also found a small H i companion 64 kpc away from the galaxy 
in projection, which adds to a previously known small stellar 
companion. No clear sign of interaction has been found so far. 

NGC 2903 is a SB d galaxy showing a symmetric strong bar 
considered typical for this class of galaxies (Laurikainen & Salo 
2002). Previous observations have shown large amounts of Ha 
emission along the bar and not only at the ends of the bar and 
nuclear region (Sheth et al. 2002). The C0(/ =1-0), Spitzer 
and GALEX data available makes this galaxy an ideal object for 
a multi wavelength study to retrieve insight in star formation his- 
tory in bars. 

A previous study by Leon et al. (2008) on the NGC 2903 
bar showed that HCN(l-O) is distributed along the bar and in the 
center. They compared the star formation rate ratio between the 
bar and the center with results from numerical simulations by 
Martin & Friedli (1997). This made them propose that the bar in 
NGC 2903 has an age between 200 and 600 Myr. 

The plan for the article is the following: in section 2 we 
present the observational data; in section 3 we analyse the gen- 
eral morphology of NGC 2903. We present our methodology to 
obtain the bar Hii regions and UV emission knots catalogues 
in section 4. The following section contains the main results re- 
garding EW//Q,, star formation rates, UV colours and ages of the 
stellar clusters. In section 6 we discuss our main results and fi- 
nally, we end with a summary and conclusions. 

2. Observations and data reduction 

2.1. Ha imaging 

The Ua data of NGC 2903 was obtained with the Wide Field 
Camera at the 2.5m Isaac Newton Telescope at the Roque de 
los Muchachos Observatory in October 2007. A summary of 
the observations is presented in Table 1. The galaxy was ob- 
served through a 95 A width narrowband filter with a total ex- 
posure time of 3600s, and an R-band filter, used for the con- 
tinuum subtraction, with an exposure time of 750s. We carried 
out the overscan subtraction, bias and flatfield correction using 
standard reduction tasks within IRAF^ . The sky subtraction was 
carried out by fitting a first order polynomial to a sky map cre- 
ated using the mean sky values derived from several positions 
on the image, free from foreground stars and galaxy emission. 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agreement 
with the National Science Foundation. 




Measured courts (x 10^) 



Fig. 1. Scaling relation from our measured H(2+[Nii] fluxes in counts 
to earlier measurements made by Mayya et al. (1994) for the same Hii 
regions. The Hii region fluxes measured by Mayya et al. (1994) also 
contain the [N ii] emission. 



The Ha-hcontinuum and R-band images were then aligned us- 
ing positions of field stars in the images and combined sepa- 
rately to produce the final images. Before continuum subtraction 
the Ha^-hcontinuum image was degraded to the resolution of the 
R-band image (~ 2.1"). 

The continuum subtraction was done as described by Relafio 
et al. (2005). We used 1 1 non-saturated foreground stars in both 
the Ha (ON-band) and in the continuum R-band (OFF-band) 
images. The fluxes of these stars yielded a mean flux ratio ON- 
band/OFF-band of 0.34+0.02. This value is then used as a start- 
ing point scaling factor to generate a set of OFF-band images 
scaled with factors ranging from 0.32 to 0.36. These images 
were afterwards subtracted from the ON-band image to produce 
images free of continuum emission. After a close inspection of 
the resulting images we adopted a value of 0.34+0.01 for the 
continuum scaling factor. The uncertainty of the continuum scal- 
ing factor produces diflTerences of up to -3% in the integrated 
fluxes of the H ii regions. 

The astrometry on the Ha image was performed using 
the USN02 catalogue coordinates for the foreground stars of 
the galaxy images, resulting in an accuracy of -0.39". The 
observations were made during non photometric conditions. 
Therefore, for the flux calibration, we used previously re- 
ported fluxes of 11 regions Hii regions located in the disk of 
NGC 2903 pubHshed by Mayya (1994). We compared the flux 
reported there with our measured flux (in ADU) using iden- 
tical apertures. We then calculated the ratio for each region 
and took the median as our final Ha flux calibration factor 
((2.2+0.4)xl0"^^erg s"^ cm"^ count"^). The error was calcu- 
lated by taking the standard deviation of the flux ratio for all 
the stars. The scaling relation is presented in Fig. 1. The 2cr sen- 
sitivity limit of the our final Ha image is 6.2 x 10"^^ erg s"^ cm"^ 
per pixel. 

The Mayya et al. (1994) Ha data from which our fluxes 
were flux calibrated were not corrected for Galactic extinction. 
Therefore, the H ii region luminosities of our catalogued H ii re- 
gions were corrected a posteriori for foreground Galactic extinc- 
tion (Schlegel et al. 1998) implying Ana =0.08 mag. 
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The bandwidth of the Ha filter also contains emission from 
the [N ii]/16548A and [N ii]A 65 84 A emission lines at their corre- 
sponding red-shifted wavelengths. A proper correction for this 
requires spectroscopic information for all the Hii regions. This 
information is not currently available. The typical [NiiJ/Hof ra- 
tio for extragalactic H ii regions of solar metallicity is approxi- 
mately 0.33 (e.g. Kennicutt & Kent 1983). Therefore, an approx- 
imate correction of 25% should be applied to the Ha fluxes and 
equivalent widths. Rather than applying this average correction 
to the fluxes of all H ii, we have opted to give the uncorrected 
Ha-F[Nii] luminosities in Table 2. 



2.2. Spitzer Infrared images 

To investigate the extinction and the star formation in the bar of 
NGC 2903 we make use of the Spitzer 24 jim observations from 
the Multiband Imaging Photometer (MIPS) instrument (Werner 
et al. 2004; Fazio et al. 2004; Rieke et al. 2004). The data was ob- 
tained in April 2004 as part of the Local Volume Legacy Survey 
(Lee at al. 2008; Kennicutt et al. 2007). MIPS reduction steps 
are described in detail by Gordon et al. (2005). The field of view 
covered by the observations is large enough to contain the whole 
galaxy with enough sky coverage allowing a good sky subtrac- 
tion. The 24 /im point-spread function (PSF) have a measured 
FWHM of -5.6" and a 2cr sensitivity limit of 0.52 x 10"^ Jy 
arcsec"^. The image was registered to the same coordinate sys- 
tem and pixel scale as the Ha image and the background was 
subtracted in the same way as for the Ha image. 

The 24 Jim image has a much larger PSF FWHM than the Ha 
image (FWHM ~2.1"). Therefore, to allow for accurate compar- 
isons, we made a copy of the Ha image degraded to the 24 jum 
PSF by convolving with a PSF kernel as described by Gordon 
et al. (2008). This copy was only used to calculate the Ha at- 
tenuation of the Hii regions (Sec. 4.3). As complementary data 
we also use the 3.6 and 8//m images, with PSF FWHM ~2" in 
both bands. These images are only used to carry out global mor- 
phological comparisons with the rest of the bands; therefore, we 
registered the images to the same coordinate system as the Ha 
image, but we did not degraded their spatial resolution. 



2.3. GALEX images 

To analyse the UV emission and its distribution in NGC 2903, 
we use GALEX mission data (Martin et al. 2005). The GALEX 
direct imaging observations are centered at 1529 A and at 2312A 
for the far-ultraviolet (FUV, 1350 - 1750A) and near-ultraviolet 
(NUV, 1750-2750A) bands respectively. NGC 2903 was ob- 
served in February 2004 as a part of the Nearby Galaxies Survey 
(NGS, Bianchi et al. 2003) with an exposure time of 1862s and 
1861s for the FUV and NUV respectively. The PSF FWHM of 
the images are ~4.2 and ~4.6" for the FUV and NUV respec- 
tively, and their respective sensitivity limits 5.25 X 10"^^ and 
3.09 X 10"^^ erg s-^cm'^A"^ per pixel (see Bianchi et al. 2003 
for further details on the NGC 2903 GALEX data). 

Both the FUV and NUV images were aligned to the same 
coordinate system and pixel scale as the Ha image and the back- 
ground was subtracted in the same way as for the Ha and Spitzer 
data. 



2.4. Other data: Ground-based optical images and 
CO (J =1-0) 

Ground based optical images from the Sloan Digital Sky Survey 
(SDSS) and a map of the CO (/=l-0) rotational transition emis- 
sion have also been used for the morphological analysis of 
NGC 2903. The optical data set comprises r\ g' and z' band 
images, which were observed as part of the SDSS data release 6 
and have a PSF FWHM of ~ 1.1". We subtracted the sky emis- 
sion from these images using the same IRAF script used for the 
Ha image (Sect. 2.1). 

For the photometric calibration, we estimated the calibration 
factor for converting the image digital counts to calibrated g\ 
/ and AB magnitudes as described on the SDSS web page 
for the data release 6^. These calibrated images were used to 
produce a ^' - colour map of NGC 2903 (Fig. 2e). 

The CO (/=l-0) emission line map of NGC 2903 has been 
used to trace the molecular gas emission and was obtained 
from the BIMA (Berkeley-Illinois-Maryl and Association) array 
and the NRAO (National Radio Astronomy Observatory)^ 12m 
single-dish telescope, as a part of the BIMA SONG key project. 
The emission line map has a maximum FWHM of 6.8". Data ac- 
quisition and reduction details can be found in Regan et al. 2001 
and Heifer etal2003. 



3. NGC 2903 general morphology 

The different appearance of galaxies in different wavelength 
ranges has already been illustrated by many authors. Galaxies 
in the UV bands seem to be patchier and later type than the same 
galaxies observed at optical and near infrared wavelengths. The 
stellar components such as bulges, bars and old stellar disks tend 
to disappear when observed in the UV range (e.g Kuchinski et 
al. 2000). NGC 2903 is a good example of this behaviour. Fig. 2 
shows NGC 2903 observed in different wavelength regimes: 
FUV, NUV, 3.6 jum 24 jum, g' - z! colour map and a CO (/=l-0) 
rotational transition emission line map. All images show over- 
laid Ha contours. 

The UV distribution is patchy and spiral-like and does not 
resemble the smooth bar-like distribution shown in the 3.6 pm 
image. The CO (/=l-0) traces the 3.6 pm. emission, while the 
Ha follows an s-shape distribution also followed by the main 
dust-lanes, as shown by the g' - z! colour image. Dust spurs 
emerge from the main dust lanes at many locations. 

The Ha leads the CO (J=l-0) emission, taking into ac- 
count that NGC 2903 rotates counter-clockwise (Hernandez et 
al. 2005), assuming trailing spiral arms. 

4. Photometry of star-forming regions 

4. t. HII region catalogue 

To study the properties of the recent star forming sites in the 
NGC 2903 bar we produced a catalogue of the H n regions lo- 
cated in the bar zone (except the nucleus). We consider the bar 
zone as the area covering the whole lense region (Kormendy 
1979.) as observed in the 3.6 pm image, see Fig. 2c. The se- 
lection criteria for considering a feature in the Ha image as an 
H II region is that the feature must have an area in pixels equal 
or larger than the image spatial resolution (i.e. area«32 pixels; 

^ http://www.sdss.org/dr6/algorithms/fluxcal.html 
^ The National Radio Astronomy Observatory is a facility of the NSF 
operated under cooperative agreement by Associated Universities, Inc. 
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(a) NUV with overlaid Ha countours. 



(b) FUV with overlaid Ha contours. 




(c) 3.6 yt/m with overlaid Ha contours. 




(d) 24 Jim with overlaid Ha contours. 





R.A. (J2000) 

(e) g'-z' colour map with overlaid Ha contours. 



(f) CO (/=l-0) with overlaid Ha contours. 



Fig. 2. Images of NGC 2903 in different spectral bands: (a) GALEX NUV image, logarithmic scale, (b) GALEX FUV image, logarithmic scale, (c) 
Spitzer 3.6yL/m, (d) Spitzer lAjim, (e) g' - z' colour-map. (f) CO (/=l-0) transition map of NGC 2903. All the colour scales are given in arbitrary 
units. All images (a) to (f) show overlaid Ha contours, the resolution of each image is given in the text. Sect 2. The Ha contour levels range 
between 4.5x10^^ erg/s/pix and 4.5x10^^ erg/s/pix (the Ha image pixel size is 0.33"/pix). 
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Fig. 3. Ha continuum-subtracted image of NGC 2903. The white open 
circles indicate the position and size of the integration apertures of the 
catalogued Hii regions, belonging to the bar region. The black open 
circles indicate the H ii regions closest to the catalogued UV complexes 
without Ha counterparts. 



with a pixel size of 0.33"/pix), having all pixels an intensity of at 
least three times the r.m.s. noise above the local background in- 
tensity level. This selection criteria implies a detectability limit 
of - 8.3 X 10^^ ergs"^. 

A total of 67 H ii regions were catalogued and their positions, 
sizes and luminosities were measured. We performed aperture 
photometry using GAIA"^. Most of the catalogued regions are 
located in the s- shape region and its neighbourhood, a region 
which is closely related to the stellar bar. A representation of the 
catalogued H ii regions can be seen in Fig. 3. 

The H II region catalogue comprises the position, size, shape 
and HcK+[Nii] luminosity of the regions, and it is presented in 
Table 2. The table contains the Hii region identification num- 
ber in column 1; the equatorial coordinates (columns 2 and 3). 
Column 4 shows the radius in arcseconds; columns 5 and 6 the 
eccentricity and position angle (measured North- West) of the in- 
tegration aperture, and column 7 shows the decimal logarithm 
of the Galactic extinction corrected Ha-h[Nii] luminosity (not 
corrected for internal dust extinction). 



4.1.1. Ha equivalent width 

The Ha EW (hereinafter EW Ha) is a measure of the strength 
of the Ha emission with respect to the continuum emission of 
the stellar ionizing cluster. Therefore, the EWna depends mainly 
on the evolutionary stage of the H ii region (decreasing as the 



^ GAIA is a derivative of the Skycat catalogue and image display 
tool, developed as part of the VLT project at ESO. Skycat and GAIA 
are free software under the terms of the GNU copyright. 



region ages), but also on the metallicity, the IME, dust content 
and ionizing photon leakage (e.g. Bresolin & Kennicutt 1999; 
Zurita & Perez 2008). 

The Ha EWs were calculated from the ratio of our measured 
Ha luminosities (not corrected for internal dust extinction) to the 
continuum emission per A. The continuum emission of the ion- 
izing cluster has been estimated from the broad band emission in 
the R-band at the position of the catalogued H ii regions. Rather 
than using our Johnson R-band image of NGC 2903, we used 
the Sloan Digital Sky Survey (SDSS) r'-band image. The lat- 
ter has a higher spatial resolution than our image, and therefore 
allowed us to better identify and measure the continuum cluster 
emission. 

When computing the Ha EW of an Hii region we are in- 
terested only in the continuum radiation emitted by the ion- 
izing cluster. However, both the flux from the ionizing cluster 
and the continuum emission from the underlying stellar popula- 
tion of the galaxy are included within the integration apertures. 
Correcting for this contamination is not straightforward and it is 
the major source of uncertainty in Hii region broad-band fluxes 
(Zurita & Perez 2008). We followed the procedure described in 
Zurita & Perez (2008) to estimate the contribution of the under- 
lying stellar population to our measured broad band fluxes by 
using two of the methods described there. The first method con- 
sist of the measurement of the local background from the median 
value within annular apertures around the H ii regions. The inner 
radius was set to the eff'ective radius of the H ii region, and the 
annulus width to approximately three times the H ii region ra- 
dius. The areas of the image corresponding to catalogued H ii 
regions had previously been masked out in the image to avoid 
contamination within the annular apertures from neighbouring 
H II regions. For some regions, the median continuum emission 
within the annulus is clearly overestimated and yields an under- 
estimation of the cluster continuum emission, which sometimes 
even results in negative fluxes. This is generally due to strong 
spatial variation of the disk continuum emission at small scales. 
Therefore, a second approach was adopted for better constrain- 
ing the local background estimates. 

The second method uses growth curves and takes advantage 
of the fact that in most H ii regions the continuum emission is 
less extended than the Ha emitting area. Therefore, the local 
background can be obtained from within the H ii region area as 
defined by its Ha emission. For each H ii region we derived ra- 
dially averaged r'-band surface-brightness profiles. For each re- 
gion we selected a radial range not contaminated by the ionizing 
cluster and determined the local background surface brightness 
from a fit to the profile in that range. 

The r'-band continuum fluxes were all corrected for the 
underlying stellar population contamination using both meth- 
ods (when possible), giving us two estimates for the r'-band 
flux of the ionizing cluster, from which we derived two val- 
ues for the EWnaHNii] • We calculated the difl'erence between 
the two EWnaHNii] estimates for each H ii region, and this was 
taken as the measured HQf-h[N ii] EW uncertainty, so it represents 
the range of values covered when using the difl'erent estimates. 
Those regions for which the two estimates yielded difl'erences 
larger than 600A, and those regions for which only one back- 
ground measurement was available, were discarded, as we con- 
sidered the results non-reliable. The resulting Ha EWs (includ- 
ing emission from the N ii lines) are presented in column 8 of 
Table 2. 
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Table 2. Hii region catalogue of the bar zone of NGC 2903. Column 1: ID number of the catalogued Hii regions; Cols. 2 and 3: equatorial coordinates; Col. 
4: Hii region radius in arcseconds; Cols. 5 and 6: Eccentricity and position angle (measured North- West) of the integration apertures; Col. 7: decimal logarithm 
of the Hof luminosity (in erg s~^); Col. 8: decimal logarithm of the Ha equivalent width; Col. 9: decimal logarithm of the 24 //m emission within the Hii region 
apertures; Col. 10 Ha attenuation see Sect. 4.3; Col. 11: estimate of the current star formation rate for each catalogued Hii region. Cols. 7,8,9 and 10 include the 
corresponding correction for Galactic extinction (Schlegel et al. 1998). The SFRs have been obtained after applying a correction of 25% for [Nii] contamination to 
the Ua luminosities and dust extinction as given in column 10. 
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4.2. UV photometry 

We performed aperture photometry on the NUV and FUV 
GALEX images to catalogue all the UV emitting knots of the bar 
region with the aim of obtaining their luminosities and colours. 
The aperture photometry has been carried out as described in 
Sect. 4.1 for the Hii region catalogue, but using circular aper- 
tures centered on the FUV peaks. The location of the apertures is 
shown in Fig. 4. The same apertures and centers have been used 
to measure fluxes on the Ha and 24 jim with the aim of deriving 
the dust attenuation (Sect. 4.3). The aperture size selection, 13.2 
arcsec, which corresponds to -570 pc at the galaxy distance, was 
determined by the spatial resolution of the 24 jim image and the 
emission peak shifts of the emitting knots at the difl'erent bands 
(see below in this section). These large apertures can cause some 
overlap between integration areas of neighbouring regions, but 
in none of the cases the contamination is significant. 

As already mentioned for the Ha photometry in Sect. 4.1, 
the background emission has to be subtracted from the mea- 
sured 24 yum, NUV and FUV fluxes within the aperture, to yield 
the stellar cluster emission alone. For estimating the background 
contamination, annuli around each aperture are difficult to define 
because of crowding. As an alternative method, we fit light pro- 
files perpendicular to the bar major axis averaged over a rectan- 
gle of width ^ 40" (much larger than our photometric apertures 
for averaging out small scale background variations, but small 
enough to ensure a valid local background estimate). The local 
background of each region is then obtained by interactive inter- 
polation of the background on the corresponding fitted profile 
(second order polynomial) at the position of the emission peaks. 
A total of 56 regions were catalogued in the bar region alone 
(excluding the galaxy center) and the positions, luminosities and 
UV colours are shown in Table 3. 

The peak emission in the UV and in the 24 jim (and Ha ) 
image is displaced a few arcseconds with respect to each other 
in most of the regions. This off'set between the emission peaks in 
UV and 24 jim or Ha amounts up to 6" and was already noticed 
by Calzetti et al. (2005) for M51. 

As already said, the 24 jum and Ha emission peaks do not al- 
ways fall within the UV apertures. Therefore, the Ha and 24 jim 
luminosities shown in Table 3 should be used with care, as they 
show the luminosities of regions centered on the UV peaks with 
the unique aim of obtaining the dust attenuation in the UV. For 
the Ha luminosities of the H ii regions we refer the reader to 
Sect. 4. 1 and Table 2. Last column of Table 3 shows the ID num- 
ber of the Hii region catalogue, matching to the UV emitting 
knot. 

A careful visual inspection of the UV, the Ha and the 24 jim 
images of NGC 2903, shows that a number of compact bright 
UV emission knots, mainly located in the bar region, in the 
Southeastern and Northwestern area from the nucleus, do not 
have significant Ha nor a 24 jim counterpart emission (see 
Fig. 2). 



4.3. Dust attenuation 

Since many years a lot of eff'ort has been put into deriving re- 
liable and easy to use dust extinction indicators to correct op- 
tical and UV observations (e.g. Kennicutt 1983; Calzetti et al. 
2000; Iglesias-Paramo et al. 2006, etc). The appearance of multi- 
wavelength observations covering from the FUV to the mid- 
infrared of a relatively large sample of nearby galaxies (e.g. 
SINGS; Kennicutt et al. 2003) has allowed to us determine reli- 
ably extinction correction using all the available bands. The Ha 



and 24 jim combination has shown to be a very useful tool to 
estimate attenuation corrected Ha and UV fluxes (Kennicutt et 
al. 2007; Kennicutt et al. 2009). We have used these results to 
obtain dust-corrected UV fluxes. We calculate the attenuation in 
the UV bands from the Ha attenuation using the laws presented 
by Calzetti (2001): 

.16, star ~ 
^0.2^,star = 

where Ax is the attenuation, in magnitudes, at the specific band. 
To obtain the Ha attenuation we used the method described by 
Kennicutt et al. (2007) : 



Ana = 2.5 log 



1 -h 



aL24u 



(LHa)obs 



where (LHa)obs refers to the observed Ha luminosity with cor- 
rection for the N ii lines, L24j2m is defined as the product vLy and 
a is the L24^^/L//q, constant scaling factor, which is fitted em- 
pirically. A value of 0.038 was adopted for a, as derived in the 
same paper for the derivation of individual region parameters, 
and a value of a = 0.02 was adopted when deriving galaxy inte- 
grated parameters. Galactic extinction was applied for the calcu- 
lation of the dust attenuation. Errors in the attenuation and UV 
magnitudes follow from two uncertainties; first, the photometric 
calibration uncertainties and second, the local background esti- 
mates. Errors are of the order of 10%. 

Because no dust tracers are detected for the UV emitting re- 
gions without Ha and 24 jum emission counterparts, we adopted 



09 32:10.0 
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Fig. 4. GALEX NUV image of NGC 2903 with the apertures used for 
photometry in the FUV, NUV, Ra and 24 jum images. See Sect. 4.2 for 
further details. 
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Fig. 5. Distribution of the Ha +[Nii] luminosity of the catalogued Hii 
regions of the bar zone of NGC 2903 from which the EWna were cal- 
culated. 

a value of zero for the attenuation in these regions. Besides the 
lack of Ha and 24 yum, no Syum (Dale et al. 2009) nor CO (/=l-0) 
emission is detected in these regions, indicating that there is very 
little gas or dust obscuring the UV emission. As an Lna upper 
limit for these regions we assumed Lna+iNin =3.8x10^^ erg s"^ 
which is the integrated Ha luminosity in a 13.2 arcsec aperture 
assuming a constant flux per pixel equal to the r.m.s. noise of the 
Ha image background. The derived Ha attenuations for the com- 
plete catalogue are presented in Table 3. We then derived UV- 
corrected magnitudes and calculated the FUV-NUV colour for 
each region. These results are also presented in Table 3: Column 
1 presents the ID of the regions; Column 2 & 3 the equatorial 
coordinates; Columns 4-7 the Ha , 24 jim FUV and NUV lu- 
minosities respectively; Column 8 the Ha attenuation; column 
9 and 10 the FUV and NUV attenuation-corrected magnitudes; 
Column 11 the FUV - NUV colour derived and Column 12 the 
H II region ID (taken from table 2) located close to each of the 
corresponding UV knot. 

5. Results 

5.1. Ha luminosity and EWna 

The log of the Ha +[NII] luminosities of the Hii regions of 
NGC 2903 range from 38.0 to 39.8 erg s"^ and their distribu- 
tion is shown in a histogram in Fig. 5. The luminosities are 
within typical values for extragalactic Hii regions in barred 
and unbarred galaxies (see e.g. Kennicutt et al. 1989; Mayya 
1994; Gonzalez-Delgado & Perez 1997; Rozas et al. 1999; 
Relano et al. 2005), but are higher than the Ha luminosities re- 
ported by Alonso-Herrero et al. (2001) for the nuclear region of 
NGC 2903. This apparent discrepancy is partly due to a smaller 
distance assumed for NGC 2903 by these authors (6.3 Mpc, ver- 
sus 8.9 Mpc assumed in this paper). More important, Alonso- 
Herrero et al. (2001) apply a semiautomated method based on a 
limiting Ha surface brightness to define the edges of an H ii re- 
gion from an HS T image, yielding catalogued Ha knots which 
may correspond to parts of a single Hii region, as defined in 
ground-based Ha imaging, and therefore have lower luminosi- 
ties and sizes. 

H II region number 36, located in the northernmost part of the 
bar (see Table 2 and/or Fig. 3), has an Ha luminosity of 6.4 x 10^^ 



erg s"\ which makes this region the brightest Hii region in the 
galaxy (outside the nuclear region of the galaxy). Its luminos- 
ity is comparable to, for instance the brightest H ii region in the 
nearby spiral galaxy M51 (Calzetti et al. 2005; Scoville et al. 
2001) and 30 Doradus in the Large Magellanic Cloud (Kennicutt 
etal. 1989). 

The EWnaHNii] of H ii regions in the bar of NGC 2903 range 
from ~ 71 to ~ 1550 A, which is typical for disk and bar Hii 
regions (e.g. Bresolin & Kennicutt 1999; von Hippel & Bothum 
1990; Cedres et al. 2005; Zurita & Perez 2008). The distribution 
of EWnaHNii] has a mean value of log EWHa+[Nii] = 2.6 ± 0.3 
(EW^.inA). 

The Ha equivalent width of H ii regions depends on several 
factors which include the initial mass function (IMF), the evolu- 
tionary status of the H ii region, its metal content, and ionizing 
photon loss due to dust extinction or leakage from the region. 
Zurita & Perez (2008) analysed the eff'ect of all these factors 
on the observed distribution of EWna for the bar H ii regions 
in NGC 1530 and we refer the reader to that paper for a more 
detailed discussion on the subject. Assuming similar physical 
conditions in the bar of NGC 2903, the evolutionary status of 
the H II region is the dominant factor aff'ecting the EW Ha dis- 
tribution. Assuming also that the H ii region is the product of a 
single burst of star formation, population synthesis models such 
as Starburst99 (Leitherer et al. 1999) can be used to estimate 
an average age for each H ii region, for a given IMF and metal- 
licity. We have therefore compared the theoretical predictions 
of Starburst99 for the evolution of a single burst of star forma- 
tion with age, with our EWna measurements. Our Starburst99 
simulations use the Geneva stellar-evolutionary tracks, and a 
metallicity of Z = 0.02 was assumed, which is the mean oxy- 
gen abundance of the gas in NGC 2903 at the range of galac- 
tocentric distances covered by the NGC 2903 bar (Pilyugin et 
al. 2004). For the IMF we took a multi power-law parametriza- 
tion (dN/dm oc m"^) with a = 1.35 for 0.1 < M < O.5M0 and 
a = 2.35 (Salpeter) for higher masses up to lOOM©. From the 
comparison of the Starburst99 output with our measurements 
(after correcting for [Nii] contamination), we estimate that the 
age of the catalogued H ii regions of NGC 2903 is within 2.7 Myr 
(log EW//^ = 3.07) and 7.1 Myr (log EW//^ = 1.73). Note that 
to derive an age from the measured equivalent width, we first 
had to correct the measured equivalent width for the [N ii] lines 
contribution. 

Figure 6 shows the distribution of the EW//^ in the galaxy 
on top of 3.6 yum intensity contours. As mentioned earlier, the 
3.6 yum emission is a strong tracer of old underlying stellar pop- 
ulations and therefore a good indicator of the stellar bar. The 
distribution of the EW//^ shows no correlation with the position 
within the bar. 

5.2. Dust attenuation 

Figure 8 presents the distribution and amplitude of the Ha atten- 
uation in the bar of NGC 2903 with the 3.6 jim intensity contours 
(tracing the stellar bar) overlaid on top of it. The attenuation was 
calculated using the ratio between the Ha emission and the 24 
yum (Kennicutt et al. 2007) as presented in Sec. 4.3. The aper- 
tures are centered on regions of strong UV emission. 
The mean attenuation value is around Ahu ~ 1 .0 mag and ranges 
from to 1.5 mag (except for regions 55 and 56, closer to the 
nucleus with an attenuation value of 2.2). There are a few UV 
regions where we assume no dust attenuation due to the lack 
of any significant Ha and 24 [im emission (see Sec. 4.3). The 
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Table 3. Summary of positions and photometry of the star-forming knots. Column 1: ID of the regions; Columns 2 & 3: right ascension and 
declination respectively; Columns 4-7: decimal logarithm of the Ha , 24 jim FUV and NUV luminosities respectively, these values have not 
been corrected for extinction; Column 8: estimated Ha attenuation in magnitudes; Columns 9 and 10: FUV and NUV attenuation-corrected 
magnitudes; Column 11: the FUV-NUV colour of the star forming region. Column 4 shows luminosities corrected for a 25% [Nii] contamination. 
For the regions where Column 8 has a value of 0.0, Columns 4 and 5 correspond to an upper limit of the corresponding luminosities. 



ID 


R.A. 

(J2000) 


Decl. 
(J2000) 


log Lna 
(erg s"^) 


log L24,um 

(erg s"^) 


log Lfuv 
(erg s"^) 


log Lnuv 
(erg s"^) 


Ana 
(mag) 


FUV 
(mag) 


NUV 
(mag) 


FUV - NUV 
(mag) 


Ha ID 


1 


09:32:11.67 


21:30:40.80 


39.15 


40.92 


40.63 


40.73 


1.27 


15.74 


15.67 


0.08 


2 


2 


09:32:12.34 


21:30:40.00 


38.93 


40.57 


40.53 


40.61 


1.06 


16.37 


16.22 


0.15 


14,15 


3 


09:32:12.25 


21:30:23.80 


38.87 


40.32 


40.80 


40.83 


0.78 


16.19 


16.03 


0.16 


4,20,21 


4 


09:32:12.75 


21:30:11.70 


<38.48 


<39.41 


40.65 


40.58 





17.95 


17.68 


0.27 


_ 


5 


09:32:12.23 


21:31:05.70 


38.94 


40.71 


40.74 


40.81 


1.29 


15.44 


15.45 


-0.01 


1 


6 


09:32:10.89 


21:30:26.10 


38.65 


40.47 


40.69 


40.83 


1.36 


15.44 


15.30 


0.14 


_ 


7 


09:32:11.78 


21:30:12.50 


38.86 


40.32 


40.90 


40.91 


0.8 


15.89 


15.81 


0.08 


2,3 


8 


09:32:11.85 


21:29:55.30 


<38.48 


<39.53 


40.77 


40.80 





17.65 


17.13 


0.52 


13 


9 


09:32:12.23 


21:30:05.70 


<38.48 


<39.75 


40.71 


40.74 





17.79 


17.26 


0.53 


51 


10 


09:32:10.18 


21:29:35.90 


<38.48 


<39.82 


40.79 


40.87 





17.59 


16.95 


0.64 


_ 


11 


09:32:10.74 


21:29:22.90 


<38.48 


<38.94 


40.41 


40.54 





18.55 


17.78 


0.77 


_ 


12 


09:32:09.78 


21:29:22.80 


<38.48 


<39.53 


40.32 


40.41 





18.77 


18.10 


0.67 


_ 


13 


09:32:07.59 


21:28:55.40 


38.97 


40.64 


40.52 


40.64 


1.11 


16.29 


16.08 


0.21 


17 


14 


09:32:07.16 


21:28:46.10 


38.75 


40.03 


40.43 


40.46 


0.6 


17.45 


17.21 


0.23 


25 


15 


09:32:08.77 


21:28:46.40 


39.17 


40.84 


40.64 


40.71 


1.1 


16.00 


15.91 


0.09 


24 


16 


09:32:07.98 


21:28:40.00 


38.88 


40.51 


40.54 


40.59 


1.06 


16.35 


16.29 


0.06 


42 


17 


09:32:07.81 


21:29:38.70 


38.90 


40.38 


40.09 


40.23 


0.83 


17.88 


17.49 


0.39 


_ 


18 


09:32:08.27 


21:29:47.70 


38.87 


40.26 


40.05 


40.14 


0.71 


18.18 


17.86 


0.33 


_ 


19 


09:32:12.03 


21:31:33.30 


38.70 


40.34 


40.54 


40.58 


1.07 


16.31 


16.31 


0.01 


26 


20 


09:32:12.12 


21:31:22.90 


38.70 


40.49 


40.64 


40.68 


1.32 


15.63 


15.72 


-0.1 


41 


21 


09:32:10.45 


21:31:14.50 


38.83 


40.70 


40.65 


40.75 


1.45 


15.36 


15.38 


-0.03 


46 


22 


09:32:11.25 


21:31:17.50 


38.78 


40.39 


40.62 


40.74 


1.02 


16.20 


15.95 


0.25 


23 


23 


09:32:09.72 


21:30:38.80 


<38.48 


<39.53 


40.19 


40.34 





19.10 


18.27 


0.82 


_ 


24 


09:32:08.78 


21:30:34.90 


<38.48 


<39.84 


40.26 


40.27 





18.91 


18.44 


0.48 


_ 


25 


09:32:08.80 


21:30:25.20 


<38.48 


<39.94 


40.17 


40.20 





19.15 


18.63 


0.52 


_ 


26 


09:32:09.63 


21:30:27.20 


<38.48 


<40.07 


40.12 


40.27 





19.28 


18.44 


0.84 


_ 


27 


09:32:08.25 


21:30:11.70 


<38.48 


<39.85 


40.09 


40.06 





19.35 


18.96 


0.39 


_ 


28 


09:32:08.45 


21:30:00.70 


<38.48 


<40.06 


40.17 


40.20 





19.16 


18.61 


0.55 


_ 


29 


09:32:11.27 


21:29:34.90 


<38.48 


<39.52 


40.37 


40.44 





18.64 


18.02 


0.62 


_ 


30 


09:32:11.47 


21:29:43.50 


<38.48 


<39.81 


40.39 


40.36 





18.59 


18.22 


0.37 


_ 


31 


09:32:07.17 


21:29:28.20 


38.83 


40.24 


39.79 


39.94 


0.74 


18.77 


18.31 


0.46 


6 


32 


09:32:09.60 


21:28:46.50 


39.16 


41.03 


40.37 


40.44 


1.45 


16.06 


16.16 


-0.09 


27 


33 


09:32:12.54 


21:30:53.60 


38.76 


40.63 


40.45 


40.59 


1.46 


15.84 


15.75 


0.09 


10 


34 


09:32:10.51 


21:31:38.80 


<38.48 


<40.07 


40.58 


40.64 





18.12 


17.54 


0.59 


28 


35 


09:32:11.05 


21:31:37.20 


<38.48 


<40.19 


40.54 


40.58 





18.22 


17.67 


0.56 


_ 


36 


09:32:10.29 


21:28:42.20 


38.62 


40.53 


40.40 


40.47 


1.53 


15.84 


15.98 


-0.14 


29 


37 


09:32:09.74 


21:31:49.20 


<38.48 


<39.48 


40.42 


40.49 





18.52 


17.89 


0.63 


_ 


38 


09:32:08.56 


21:31:34.30 


38.76 


40.44 


40.61 


40.65 


1.12 


16.05 


16.06 


-0.01 


30,52 


39 


09:32:07.94 


21:31:36.10 


38.82 


40.01 


40.59 


40.61 


0.5 


17.20 


16.97 


0.23 


37 


40 


09:32:07.72 


21:31:25.80 


<38.48 


<39.31 


40.54 


40.60 





18.22 


17.63 


0.59 


_ 


41 


09:32:10.56 


21:28:59.90 


38.92 


40.38 


40.67 


40.72 


0.81 


16.48 


16.28 


0.19 


_ 


42 


09:32:11.30 


21:28:56.40 


38.90 


40.76 


40.64 


40.66 


1.42 


15.46 


15.64 


-0.18 


31 


43 


09:32:11.61 


21:29:06.70 


38.75 


40.52 


40.47 


40.51 


1.27 


16.15 


16.20 


-0.05 


45 


44 


09:32:10.89 


21:28:30.90 


<38.48 


<40.03 


40.52 


40.52 





18.28 


17.82 


0.47 


_ 


45 


09:32:11.72 


21:28:31.70 


38.74 


40.16 


40.39 


40.36 


0.76 


17.26 


17.24 


0.01 


32 


46 


09:32:10.59 


21:28:21.00 


<38.48 


<39.53 


40.55 


40.57 





18.20 


17.71 


0.49 




47 


09:32:09.78 


21:28:14.80 


<38.48 


<39.74 


40.64 


40.60 





17.97 


17.61 


0.36 




48 


09:32:12.68 


21:28:37.90 


39.01 


40.57 


40.54 


40.53 


0.94 


16.55 


16.59 


-0.04 


33 


49 


09:32:13.30 


21:28:46.40 


39.07 


40.46 


41.10 


41.04 


0.72 


15.54 


15.60 


-0.06 


34,43 


50 


09:32:12.83 


21:29:38.40 


38.48 


39.71 


40.66 


40.59 





17.92 


17.64 


0.28 




51 


09:32:13.41 


21:29:49.60 


<38.48 


<39.63 


40.40 


40.33 





18.57 


18.29 


0.28 




52 


09:32:06.38 


21:31:15.80 


38.65 


39.77 


40.42 


40.44 


0.43 


17.75 


17.47 


0.28 




53 


09:32:13.65 


21:28:58.00 


<38.48 


<39.77 


40.47 


40.49 





18.39 


17.90 


0.49 




54 


09:32:09.14 


21:28:11.60 


<38.48 


<39.98 


40.48 


40.47 





18.38 


17.95 


0.43 




55 


09:32:08.05 


21:28:15.20 


<38.48 


<40.14 


40.30 


40.40 





18.82 


18.11 


0.71 


35 


56 


09:32:09.67 


21:31:08.80 


39.42 


41.32 


40.25 


40.44 


1.5 


16.26 


16.08 


0.18 


36 



higher attenuation values are located in the bar and at the begin- and spiral arms in NGC 2903 are gas rich, whereas the UV emit- 
ning of the spiral pattern (see Table 3). This shows that the bar ting regions to the Northwest and Southeast of the galaxy nu- 
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cleus, those regions with zero attenuation, have low gas and dust 
content; nevertheless, this region is populated by young stellar 
clusters, see Sec. 5.4. Further discussion on the origin of these 
regions is presented in Sec. 6. 

5.3. Star formation rate 

The locations of the most active star formation sites within a 
bar are determined by the favored dynamical behavior of gas 
and dust. Studying the distribution of the loci of current star- 
formation within bars can therefore put constraints on the prop- 
erties of a particular bar potential. Gas does not stay in typi- 
cal bar intersecting orbits and tends to accumulate and shock 
in certain regions of the bar. Dust lanes can be used as tracers of 
shocks in the gas flow (Athanassoula 1992). From Fig. 2 one can 
already see that the dust lanes (as seen from the g' - z' colour 
map) follow the CO (J=l-0) distribution, with the Ha leading 
the CO (J=l-0) as it has also been previously observed in other 
barred spirals (e.g. Sheth et al. 2002). 

The Ha emission in the bar of NGC 2903 follows an s-shape 
covering the whole stellar bar and it is dominated by classical 
localized H ii regions (Fig. 3). 

In ionization-bounded H ii regions, the reddening corrected 
Hof luminosity scales directly with the total Lyman continuum 
flux emitted by the ionizing embedded star cluster. Therefore the 
Ha emission is a good tracer of current star formation. We have 
used the following calibration relating the Ha luminosity to the 
rate of star formation (Kennicutt 1998) 



EW^r^ Distribution 



SV^Ha (Mq year"^) = 7.9 x 10 



^-42 



L//a(ergs ^), 



for calculating the current SFR associated to each catalogued 
Hii region (Sect. 4.1). The results are given in column 11 of 
Table 2 and have been obtained from the observed Ha luminosi- 
ties (column 4 of Table 2), and have been corrected for [Nii] 
emission lines contamination (see Sect. 2.1) and dust extinction 
inside the regions. The latter have been estimated as already de- 
scribed in Sect. 4.3 for the UV emitting knots, comparing the Ha 
and 24 jim luminosities within the Hii region apertures. Given 
the difl'erent spatial resolution of the two images, we applied an 
aperture correction factor (Table 1 of Engelbracht et al. 2007). 
The SFRs obtained for our catalogued H ii regions ranges from 
~ 1 . 1 to ~ 126 X 10"^Mo year"^ We have excluded regions 55 and 
56 in this calculation due to the likely contamination from the 
nuclear region emission. 

The total SFR implied by the Ha emission of the bar zone 
of NGC 2903 (excluding the nucleus) is ~ 0.9 ± 0.2 MQyr-^ 
(Fig. 7 shows the integration area). This value implies that the 
NGC 2903 bar is forming stars at a rate comparable to other 
normal (non-starburst) spirals. The nucleus of NGC 2903 is al- 
most as bright in Ha as the whole bar, implying similar SFRs 
(Alonso-Herrero et al. 2001; Leon et al. 2008). 

Other SFR calibrations make use of the UV emission. We 
have derived the integrated SFR, within the same region as for 
the calculation of the SFR//^, using the UV emission as SF in- 
dicator. The SFR in the UV has been computed following the 
recipes given in Salim et al. (2007) 

SFRuv (Moyear"^) = 1.08 x IO'^^Lfuv, 

with Lpuv in ergs s"^Hz"^ The star formation rates calculated 
in this way give a SFR/7f/y = 0.44 ± 0.06 M© yr"^ for the bar 
region. The SFR calculated from both bands (cf. SFR^/y = 
0.44 ± 0.06 Mo yr-i and SFRna = 0.9 ± 0.2 M© yr'^) and the 
fluxes from which these SFR have been obtained, are compatible 



^ +21°30'.O 




Fig. 6. Ua +[Nii] EW distribution in the bar region of NGC 2903. The 
size of the circles represents the H ii region projected area. The intensity 
contours trace the 3.6 /dm emission. 



SFR Distribution 




09:32:20.0 15.0 09:32:10.0 05,0 09:32:00.0 

R.A. {J2000) 

Fig. 7. SFR distribution in the bar region of NGC 2903. The size of the 
circles represents the Hii region projected area. The SFR is given in 
Moyr"^ The SFR has been estimated from the measured Ha luminos- 
ity of the Hii regions (see Sect. 5.3). The intensity contours trace the 
3.6 jum emission. The solid black line indicates the integration area for 
the calculation of the total SFR in the bar region (see Sect. 5.3) 



with the results given by other works based on a large number 
of galaxies (i.e. Lee et al 2009; Salim et al 2007). The SFRuv is 
a tracer for recent star formation, averaged over some hundreds 
of Myrs (Calzetti et al 2005; Kong et al. 2004) while the SFR//^ 
probes the current star formation. The SFRs derived in this pa- 
per are compatible with a number of star formation histories, as 
derived from the models (Starburst99). It is, however, interest- 
ing to date the age of the UV complexes with no Ha counterpart 
to shed some light into the origin of these structures. This will 
be investigated in the following section. 
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Fig. 8. Distribution of the Ha attenuation in NGC 2903. The circle size 
represents the integration apertures used for the UV catalogue. The in- 
tensity contours trace the 3.6 yum emission. 
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Fig. 9. Distribution of the UV colours in the NGC 2903 bar and its sur- 
rounding regions. The intensity contours trace the 3.6 /um emission. A 
concentration of redder regions is clearly visible to the Northeast and 
Southwest of the nucleus. These regions correspond to the UV regions 
with no significant Ka counterpart. 



5.4. UV colours and cluster ages 

There is now compiling observational evidence showing that the 
FUV - NUV colours are very sensitive to the age of the young 
stellar populations due the rapid evolution of the most massive 
stars (i.e. Bianchi et al. 2006; Calzetti et al. 2005). To estimate 
the ages of our catalogued UV regions we use the FUV - NUV 
colours. FUV - NUV colours were calculated using fluxes cor- 
rected for dust attenuation (Sect. 4.3). Figure 10 shows the range 
of UV colours in the bar of NGC 2903. The colours have a mean 
value of 0.3 mag with a standard deviation of 0.27 mag. This 
range of values is typical for spiral galaxies (Thilker et al. 2005; 
Koribalski & Lopez-Sanchez 2009). Figure 9 shows the spatial 
distribution of the derived FUV - NUV regions. One of the most 
striking facts on the colour distribution are the redder regions 
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Fig. 10. FUV - NUV colour distribution for the catalogued UV regions 
in the bar zone of NGC 2903. 
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Fig. 11. Time evolution of the UV colour for a single starburst region. 
The model assumes Padova stellar tracks with Z=0.02 (solar) and an 
IMF constructed by a multi power-law parametrization with or = 1.3 
for 0.1 < M < 0.5Mq and a = 2.35 for higher masses up to lOOM©. 



located to the Northwest and Southeast of the galaxy nucleus 
and at the beginning of the spiral arm. These regions correspond 
mostly to those with the lowest attenuation. The regions located 
on the bar consist mostly of blue UV emitting regions. 

As well as age-dependent, the FUV - NUV colours are also 
strong tracers of metallicity; however, we can assume an ap- 
proximately constant metallicity in the bar region of NGC 2903 
(Pilyugin et al. 2004). With this assumption the FUV - NUV 
colour is, as we already mentioned, a good tracer of age. To com- 
pare our results with the theoretical models in order to derive the 
ages of the clusters, we model the evolution of a single star for- 
mation burst using the Starburst99 population synthesis model 
(Leitherer et al. 1999). As input values we used the Padova AGB 
stellar tracks with a metallicity of Z - 0.02 (solar). For the IMF 
we took a multi power-law parametrization with a = 1.3 for 
0.1 < M < 0.5Mo and a = 2.35 for higher masses up to lOOM©. 
We then derived FUV - NUV colours out of the simulated spec- 
trum at multiple time-steps from 1 Myr up to 1 Gyr (see Fig. 11). 
In order to estimate the UV region ages, we interpolated the 
modelled FUV - NUV colour versus age dependence at the ob- 
served UV colours. The colours of the bluer regions range from 
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-0.18 to 0.15 implying an age of 3 to 10 Myr. The redder regions 
have colours ranging from 0.4 to 0.85 which imply ages ranging 
from ~ 150 to 320 Myr. The latter group of colours corresponds 
to regions with no significant Ha nor 24 jim emission. The FUV - 
A^f/y colour remains constant around 0.2 implying ages between 
10 and a little less than 100 Myr. 

6. Discussion 

NGC 2903 is a clear example of a galaxy with very differ- 
ent morphologies as shown by different broad-band images. 
Although its bar is classified in the NIR with a bar strength of 3 
(Laurikainen & Salo 2002), corresponding to a relatively strong 
bar, and it is considered a classical bar, its morphology at shorter 
wavelengths is more similar to that of a patchy spiral disk. This 
clearly shows that to understand the SF history in bars it is nec- 
essary to have a panchromatic view. 

To investigate the properties of SF in the NGC 2903 bar, 
we have at our disposal information on the distribution of gas 
in two different phases, namely the ionised gas and the colder 
CO {J =1-0). As for the stellar component, we are able to trace 
young stars up to 1 Gyr with the UV information, the location 
of massive stars can be traced from the location of H ii regions, 
whereas intermediate and older stellar populations can be probed 
through their distribution in 3.6 [im and optical bands. 

We observe that the CO (/ =1-0) distribution follows closely 
the emission in the 3.6 yum band. There is a misalignment be- 
tween the major axis of the bar as derived from ellipse fitting 
of optical images (Sheth et al. 2002) and the 3.6 [im major axis 
light distribution. So, both the CO (/ =1-0) and the 3.6//m would 
be on the leading side of the optical photometric axis of the bar. 
The Ha emission not only leads both the CO (J=l-0) and the 3.6 
lim emission but has a slightly different distribution, following a 
spiral pattern not so clearly seen in the CO (J=l-0) and 3.6 jjm 
emission. Also, the Ha emission leads the dust lanes as traced 
by optical colour maps. Numerical models predict the gaseous 
component (both ionised and neutral gas) to lead the stellar bar 
(Martin & Friedli 1997). 

This study started as a follow-up of the analysis carried out 
on the bar of NGC 1530 (Zurita & Perez 2008) for understand- 
ing the physical processes determining star formation in bars, 
where it was clearly shown the importance of the bar dynamics. 
CO (J=l-0) and dust is concentrated on regions where shocks 
and gas accumulation is favoured by the bar potential. The dis- 
tribution of the dust in the bar of NGC 1530 is much narrower 
than the distribution in NGC 2903. In NGC 1530, the dust lanes 
are straight and very well delineated as opposed to the broken 
feathery morphology of the dust lanes in NGC 2903 (Fig. 2e). 
This fact might be indicating that no strong shocks are present 
in the bar of NGC 2903, contrary to the large velocity gradi- 
ents perpendicular to the bar in NGC 1530 (Zurita et al. 2004), 
suggesting that shocks are driving the observed morphology in 
NGC 1530. 

Streaming motions along the bar of NGC 2903 were detected 
using HCN(l-O) as the kinematic tracer by Leon et al. (2008). 
An undergoing analysis of the velocity gradients in the ionised 
gas of NGC 2903, using the data published by Hernandez et al. 
(2005), will clarify whether strong velocity gradients are present 
or not in the bar of NGC 2903. Therefore, although the NIR 
emission suggest that NGC 2903 hosts a strong bar, the broken 
and wide structure of the dust lanes might be an evidence of the 
opposite. 

The presence of H ii regions in both the trailing and leading 
sides of the dust lanes in NGC 1530 and the evidence of an age 



gradient between both sets of H ii regions, gave support to the hy- 
pothesis (Sheth et al. 2002) that massive stars form in the trailing 
side of the bar dust lane and move to the leading side. There are 
no H II regions on the trailing side of the bar of NGC 2903, and 
no trend in H ii region age is found across the bar. This may in- 
dicate that the explanation for where star formation is favoured, 
is not as straightforward as in NGC 1530, or at least it is not as 
clearly driven by the strong bar dynamics. 

The complexity of the SF properties in the bar of NGC 2903 
is also shown by the presence of compact UV sources without 
significant 24 iim and Ha emission. These regions are located 
symmetric to the center of the galaxy and nearly perpendicular 
to the bar as well as in an spiral pattern. All of these regions have 
FUV-NUV colours compatible with ages of a couple of hundred 
Myr. The lack of Ha emission suggests that these locations are 
currently not favoured by the bar potential to form stars. These 
stars could have formed somewhere else in the bar and migrated 
afterwards to the current position. But then a few questions arise: 
where have these stars formed?; are the stars currently forming 
likely to populate similar locations? or, on the contrary, were 
they formed during a more general SF burst possibly linked to 
the bar formation, or a merger? An undergoing analysis on the 
H I and Ha velocity fields will shed some light into the origin of 
the UV complexes. 

The current SFR in the bar, as calculated from Ha, shows 
values similar to other disk galaxies although it has considerable 
high amount of dense cold gas, as traced by the HCN(l-O) (Leon 
et al. 2008). The recent SFR, as calculated from UV emission, is 
compatible with that of the current SFR. 

It is worth noticing the offsets of a few arcseconds between 
the emission peaks in the 24 jjm and the GALEX images. 24 //m 
is a good tracer of current star formation while UV emission is 
a 'time averaged' recent star formation indicator (e.g Calzetti et 
al 2005). The offsets found here (also seen in M51 by Calzetti et 
al 2005) show an evolutionary link between the sites where cur- 
rent star formation is occurring and the position of the UV coun- 
terpart. Estimations of the rotation curve and the pattern speed 
would be necessary to determine if the dynamical time-scales 
are compatible with the aging of the regions. 



7. Summary and conclusions 

We have performed a detailed multi- wavelength study from UV 
to sub-millimeter observations on the NGC 2903 bar and its sur- 
rounding regions. We mapped and catalogued the H ii regions of 
the bar and measured their Ha equivalent widths. Furthermore 
we have obtained a catalogue of the UV emitting regions, and 
measured their peak location with respect to Ha . The extinction 
has been estimated using the Ha and 24 fim emission. SFRs us- 
ing both Ha and UV indicators have been calculated. We have 
estimated the age of the regions, using the EWna and the FUV- 
NUV colour together with stellar population synthesis models. 
Our main results are: 

- NGC 2903 is a morphologically complex galaxy. The near- 
infrared as well as the CO (7=1-0) band show a clear barred 
structure whereas the Ha and UV maps show a patchy spiral 
like structure. 

- There are clear spiral like UV complexes with no signifi- 
cant Ha, 24 ^m and CO (J=l-0) counterpart emission. These 
complexes are located Northwest and Southeast of the bar 
within the inner 1 arc minute radius (corresponding to ~2.5 
kpc). 
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- The Ha emission along the bar, leads the CO. The 3.6 //m 
and CO (J=l-0) emission trace each other, both leading the 
major axis of the optical light distribution. 

- The Ha luminosities and EWna of the bar H ii regions are 
within typical ranges for H ii regions in bar and unbarred spi- 
rals. The spatial distribution of the Ha EWs does not corre- 
late with any morphological feature of the bar. 

- The average dust attenuation in the bar area of NGC 2903 is 
AHa= 1-06, and ranges from to 1.5 mag. 

- The FUV-NUV colour distribution is distributed in two re- 
gions, the bluer regions range from -0. 1 8 to 0. 15 implying an 
age of ~3 to 10 Myr. The redder regions have colours rang- 
ing from 0.4 to 0.85 which imply ages ranging from 150 to 
400 Myr. The latter correspond to regions with no significant 
Ha nor 24 jim emission. 

- The SFRs of the bar region derived from Ha and from the 
UV emission are 0.9 ± 0.2 MQyr~^ and 0.4 ± O.lMojr"^ 
respectively. 

All these results suggest that an agent triggered a SF burst a 
few hundred Myrs ago. Interestingly, we see some stellar clus- 
ters unrelated to the current SF locations, symmetrically located 
nearly perpendicular to the bar (in the inner ~2.5 kpc). The ori- 
gin of these regions might be related to this SF burst. 

In a following paper we will analyse the gas kinematics of 
the galaxy, with the aim of shedding some light on the origin 
(merger vs. secular evolution) of these findings. 
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